clear

global b;
global I;
global Kd;
KArr = [0:0.01:0.1];
%K = KArr(1);
Kd = 0.011;
b = 0.02;
I = 0.004;

%for i=1:length(KArr)
    %Kd=KArr(i);
    [T Z]=rk4sgral(0,1000,[0,0],10000,'thetapuntoe');
    plot(T', Z(:,1), 'r-');
    %, T', Z(:,2), 'b^'
    hold on;
%end

hold off;

%res(:,1), thetaR(res(:,1)), 'b-');